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Abstract. A quantitative analysis is conducted on the impacts of experimental 
5 ' imperfections in the input state, the detector properties, and their interactions on 

r^ , photon-subtracted squeezed vacuum states in terms of a quantum non-Gaussian 

^^' character witness and Wigner function. Limitations of the non-classicality and 

quantum non-Gaussian characteristic of Schrodinger kitten states are identified and 
addressed. The detrimental effects of a photon-number detector on the generation 
of odd Schrodinger kitten state at near-infrared wavelengths (^ 860 nm) and 
telecommunication wavelengths (^ 1550 nm) are presented and analysed. This 
^^ ' analysis demonstrates that the high dark count probability of telecommunication- 

^D ■ wavelength photon-number detectors significantly undermines the negativity of the 

'!::j- I Wigner function in Schrodinger kitten state generation experiments. For a one-photon- 

^P ' subtracted squeezed vacuum state at ~ 1550 nm, an APD-based photon-number- 

resolving detector provides no significant advantage over a non-photon-number- 
resolving detector when imperfections, such as dark count probability and inefficiency, 
are taken into account. 



1. Introduction: 

Non-Gaussian operations on states have attracted intense interest in quantum 
continuous variable (QCV) information processing as they provide significant advantages 
for universal quantum computing [1], quantum teleportation [2], entanglement 
distillation [3, 4], high- precision measurement [5], and proposed loophole- free tests of 
Bells inequalities [6]. Nowadays, two main categories of non- Gaussian and non-classical 
optical quantum states with negative-valued Wigner functions, such as Fock states 
[7, 8, 9] and photon-subtracted squeezed states [10, 11, 12, 13], have been experimentally 
generated based on parametric down- conversion (PDC) in nonlinear crystals followed by 
photon number detection. In contrast to non-degenerate PDC for Fock state generation, 
degenerate PDC is used to generate photon-subtracted squeezed states. A small fraction 
of the squeezed vacuum beam is tapped off via a beam splitter and guided into a 
photon-number detector. The tapped-off light is used as a trigger to condition the 
remaining beam into a photon-subtracted squeezed vacuum state [11]. The projected 
state is referred to as an optical Schrodinger kitten as it closely approximates an optical 
Schrodinger cat state with a small amplitude [10, 11]. 

Photon-subtracted squeezed vacuum states with negative Wigner functions have 
been successfully demonstrated using Ti: Sapphire lasers at wavelengths around 860 
nm and nonlinear crystals such as potassium niobate (KNb03) and periodically- 
poled KTiOP04 (PPKTP) [10, 11, 12, 13]. However, to our knowledge, a negative- 
valued Wigner function for Schrodinger kitten state at telecommunication wavelengths 
~ 1550 nm is yet to be experimentally demonstrated [14]. We hypothesize that 
the principal difference in Schrodinger kitten states generation between 860 nm 
and 1550 nm lies in the performance of photon- number detectors used in state 
preparation, which may undermine the negativity of the Wigner functions. Non- 
Gaussian states at telecommunication wavelengths are indispensable for secure optical 
quantum telecommunication due to their low loss in optical fibres. Therefore, it is 
imperative to develop a model to analyse the properties of Schrodinger kitten states, 
and improve the experimental design based on the theory. 

Historically, negativity in the Wigner function has been the standard criterion to 
identify whether a state generated from an experiment is non-classical [15, 16]. However, 
the negative Wigner function of a quantum state generated from an experiment easily 
degrades and becomes positive since the quantum state is fragile to any loss before 
it arrives at the verifying detector. Therefore, it is not suitable to solely rely on 
this criterion to characterise the non-classicality. To resolve this confiict, Jezek et 
al [17, 18] proposed a quantum non-Gaussian character witness to verify states with 
positive Wigner functions that cannot be prepared by merely using Gaussian states and 
operations. 

Dakna et al proposed the concept of conditional measurement based on a lossless 
beam splitter to generate a Schrodinger kitten state by subtracting photons from a 
squeezed vacuum state. They developed a model by taking the squeezing level, beam 



splitter transmission, photon-number detector inefficiency, and non-photon-number- 
resolving ability into account [19, 20]. The impacts of these factors on the Wigner 
function and non-classicality of a state were investigated by Olivares et al [21]. Kim et 
al analysed the necessary conditions to obtain a negative Wigner function for a realistic 
case including the input as a mixed state, threshold detection, inefficient homodyne 
detection, and mode purity in the subtraction path [22] . The dark count influence of an 
on-off photon-number detector on the Wigner function of a Schrodinger kitten state was 
considered by Suzuki et al [23] , and all of these imperfections were incorporated into a 
model developed by Gerrits et al [13]. However, a quantitative analysis has not been 
conducted of the impacts all these experimental imperfections have on the negativity of 
the Wigner function and on the quantum non-Gaussian character witness of Schrodinger 
kittens. Particularly, a thorough comparison between kitten state generation at ~ 860 
nm and generation at ~ 1550 nm is yet to be explored. 

In this paper, we quantitatively analyse the impacts of experimental parameters 
involving the impurity of the input state, inefficiency of the photon-number detector, 
dark count probability, non-photon-number-resolving ability, mode purity, and the 
homodyne detector inefficiency on Schrodinger kittens generation in terms of the 
quantum non-Gaussian character witness proposed in references [17, 18] and the Wigner 
function value at the origin, W(0,0). The comparisons between Schrodinger kittens at 
near-infrared wavelengths (~ 860 nm) and telecommunication wavelengths (~ 1550 
nm) are discussed, and principal limitations of state generation at telecommunication 
wavelengths are identified. 

The paper is organized as follows. In Section 2, a Schrodinger kitten state 
model is derived covering all possible experimental imperfections based on conditional 
measurement. In Section 3, the quantum non-Gaussian character witness introduced 
in references [17, 18] is described. In Section 4, we present the physical mechanism for 
each experimental imperfection and quantitatively analyse the effects of these impacts 
on Schrodinger kitten state features. Concluding remarks are given in Section 5. 

2. Schrodinger kitten state generation based on conditional measurement 

A schematic diagram of a theoretical model for Schrodinger kitten state generation with 
experimental imperfections is shown in figure 1 [13]. The model includes three parts: 
input state, photon subtraction, and state characterisation. The photon subtraction 
unit is composed of a 'magic' reflector that is arbitrarily tuneable via a half-wave 
plate and a polarization beam splitter [11, 13]. Ideally, when an even (odd) number 
of photons are subtracted from a pure squeezed vacuum state, an even (odd) kitten 
with a negative Wigner function can be obtained. However, numerous factors can 
undermine the ability of such an experiment to produce a Wigner function with 
negativity. These factors include optical elements related to the experiment, such as 
the impurity of the input squeezed vacuum state, mode impurity before the photon- 
number detector, and inefficiency of the homodyne detector used to characterise the 



quantum state. Imperfections in the photon-number detector can further degrade the 
prepared state. These imperfections involve a high dark count probabihty, low quantum 
efficiency, and the non-photon-number-resolving ability of some detectors. Therefore, 
a quantitative analysis of all these imperfections can shed light on the practical 
generation of Schrodinger kitten states, particularly with regards to experiments at 
telecommunication wavelengths. 
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Figure 1. Theoretical model for Schrodinger kitten state generation with experimental 
imperfections. Pin2i — Pti, fi'- reflectivity of BSl, ii: tansmission of BSl, r2: 
reflectivity of PBS, Pmii'- input state density matrix of BSl, pti'. transmission state 
density matrix of BSl, pri'- reflected state density matrix of BSl, Pin2i' input state 
density matrix of PBS, Pout2i'- output state density matrix of PBS, rjHD'- homodyne 
detection efficiency, BSl: beam splitter for modelling the impurity of the squeezed 
vacuum state, HWP: half-wave plate, PBS: polarization beam splitter, PND: photon- 
number detector. Vac: Vacuum state, BS2: beam splitter for modelling the inefficiency 
of the homodyne detector. The HWP and the PBS comprise a 'magic' reflector. The 
HWP and the PBS comprise a 'magic' reflector. 



2.1. Impure input state 

Parametric down-conversion based on a second-order nonlinear process in nonlinear 
materials is an effective approach to generate squeezed vacuum states. Ideally, a 
squeezed vacuum state consists of a photon number distribution with only even photon 
numbers. Impurity in the state contaminates the photon number distribution with odd 
photon number probabilities. This impurity can be equivalent to loss in a pure squeezed 
vacuum state, and can be described as a pure squeezed vacuum state followed by a beam 
splitter as shown in figure 1. 

A pure squeezed vacuum state with a squeezing angle of zero degree can be 
expressed as [24]: 
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where U is the unitary operator of the beam sphtter. 

The transmitted density matrix after BSl (i.e. the impure squeezed vacuum state) 
can be obtained by tracing the output density matrix over ri [25], 



where 



{2n)\{2h)\ 



Ph = Trr^ [Poutl] 

oo 2min{n,b) 

^.5o S \{2n-k)\{2h-k)\ 



Poutll — |^OUtll)(^OMtll| 



a2na2br1ei+'-'' 



k\ 



2n-k){2b-k\,{5) 



(6) 



2.2. Conditional measurement based on a lossless beam splitter 

According to the conditional beam sphtter operator in reference [20], we have 
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where {| {| denotes the magnitude of a state vector and the non-unitary beam sphtter 
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where pin2i = pti as shown in figure 1 and expressed in (5), and Tr[ ] denotes the trace 
of a matrix. 



2.3. Schrodinger kitten state prepared with an ideal photon-number-resolving detector 

In the case of an ideal photon-number-resolving detector (PNRD) (i.e. no dark counts 
and the quantum efficiency is 100%), we can obtain the projected state density matrix 
by substituting (5) into (9), 



where 
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However, an ideal PNRD is unavailable in practical experiments. Avalanche 
photodiodes (APDs) are usually used as photon-number detectors, where Si-APDs 
and InGaAs-APDs are used to detect near-infrared wavelengths (~ 860 nm) and 
telecommunication wavelengths (~ 1550 nm), respectively. Therefore, it is imperative 
to consider all possible imperfections of the photon-number detector, including the dark 
count probability, quantum efficiency, and the non-photon-number-resolving ability, and 
implement a quantitative analysis on the impact of all these experimental imperfections 
on the resultant quantum state. 

2.4. Schrodinger kitten state prepared with an imperfect photon-number detector 

2.4.1. Dark counts probability and quantum efficiency On the one hand, the existence 

of dark counts causes 'false' clicks even if a photon is not actually subtracted. On 

the other hand, some actual clicks are missed due to the inefficiency of the detector. 

Therefore, m-click events may originate from m-1, m-2,...0 or m-hl, m-l-2,... actual 

photons being subtracted. Consequently, the conditional state is a statistical mixture, 

which can be expressed as [13, 19]: 

00 
PiMPNBDim) = Y^ Qik\m)pouti{k), (12) 

fe=0 

where Q{k\m) is defined as the conditioned probability, with which m photons 
would have been subtracted, given that k photons are actually detected by the imperfect 
detector. According to the Bayes rule, we can obtain the conditional probability, 
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where Pdc and tjapd are the dark count probability and quantum efficiency of the 
detector, respectively [26]. 

2.4-2. Non-photon-number-resolving ability Most photon-number detectors used in 
experiments so far are on-off or non-photon-number-resolving detectors (NPNRDs) 
without the capability to distinguish the specific number of detected photons. Different 
from PNRD, NPNRD accepts k clicks even though the actual number of clicks can be 
larger than k. Thus we have [13] 

piMNPNRDim) = }^ ^oo n(i.\ • ^^^^ 

k=m l^k=m Vlft-J 

2. 5. Mode purity of subtracted photons 

Mode purity, s , is defined as the probability that the detected photons by the photon- 
number detector are mode matched to the local oscillator (LO) used in the kitten state 
characterisation via homodyne detection. As it is quite difficult to obtain a perfect 
mode purity, the detected density matrix of a projected state would be a mixed state 
consisting of the actual projected state and the unprojected state (i.e. the input state 
with loss). Therefore, we have [22, 28] 

Pdetect = S * Pprojected + (1 — S J * Pinputwithloss- \^') 

2.6. Schrodinger kitten state characterised by an inefficient homodyne detector 

Homodyne detection is a typical approach used to characterise the projected state. The 
homodyne detector efficiency is calculated as [11, 12, 13] 

Vhd = Vqe *Vt* C^ (18) 

where tjqe is the quantum efficiency of the two photodiodes in the homodyne detector, rjt 
is the transmission coefficient from the 'magic' reflector to the homodyne detector, and 
C is the visibility of interference fringes between the signal and LO, denoting the degree 
of mode matching. Therefore, the total efficiency, tjhd, quantifies various categories 
of loss. As show in figure 1, the homodyne detection inefficiency can be simulated by 
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a lossless beam splitter before a perfect homodyne detector, and the density matrix 
measured with an inefficient homodyne detector is given by [29], 

oo 
{l\Pdetect{VHD)\n) = ^ Bi+kAv)Bn+k,n{v) {^ + k\Pdetect\n + k) (19) 

fe=0 

in terms of the initial field density matrix pdetect, where 



Bi^kAv) = f-^if^^Hnil - Vhd)'. (20) 

3. Character witness to identify the quantum non-Gaussian state and 
non-classical state 

Up until very recently, negativity in the Wigner function has been widely used as the 
criterion to identify the non-classicality of a state [15]. However, for some non-classical 
quantum states, such as squeezed vacuum states, this criterion does not work because 
they possess positive Wigner functions. In addition, some heralded quantum states 
have positive Wigner functions that could not be prepared from Gaussian states and 
linear optical devices. Therefore, Jezek et al [17, 18] proposed a non-classical and a 
quantum non-Gaussian witness. States beyond a convex set of stochastic mixture of 
coherent states are defined as non-classical states. Similarly, quantum non-Gaussian 
states are referred to as states beyond a convex set of stochastic mixture of Gaussian 
states [17, 18]. 

The quantum non-Gaussian character witness is based on Fock state basis, and 
is introduced as a linear combination of zero photon probability, po? and one photon 
probability, pi, in the Fock state basis density matrix [17, 18], 



W{a) = apQ+pi (21) 



where 



p—e^sinhr 



Po = 7—, (22) 

coshr 

4r 1 —e^sinhr 

P^ = A /.3 ' (23) 

a G [0, 1] is a dimensionless number and r G [0, oo) is the squeezing parameter. A 
quantum Gaussian boundary, Wa{a), is defined as the maximum value of W(a) over 
a and r. The quantum Gaussian character witness value is defined as W(a) - Wda). 
If this witness value is larger than 0, then the state is a quantum non-Gaussian state. 
For quantum states related to squeezed states, such as squeezed single photon states 
or Schrodinger kitten states from photon-subtracted squeezed states, the quantum non- 
Gaussian character witness is generalized by an ant i- squeezing operation [18], 

W{a,s) = apo{s)+pi{s) (24) 
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where 

p^is) = {n\S-^{s)pSis)\n) (25) 

S~^ and S are the anti-squeezing and squeezing operators, respectively, and p 
corresponds to the density matrix of the state in the Fock state basis. The quantum 
non-Gaussian character witness value is defined as W{a,s) — WG{a,) in the following 
sections. 

Equivalently, a classical boundary is defined as the maximum value of 

WM = apQ + pi (26) 

over a, where 

Po = e-^ (27) 

pi = ne"", (28) 

and n G [0, oo] is the mean photon number. 

Therefore, it is easy to identify the quantum non-Gaussian or non-classical 
characteristic of a state via its density matrix. 

4. Quantitative analysis of experimental imperfections impact on 
Schrodinger kitten state generation 

4.I. Estimation of pure input state level and input state impurity 

Generally, squeezing value is referred to as the noise variance of a squeezed state, which 
is related to the degree of squeezing ,^ in (1) by 

Vsqz = cosh{2^)-sinh{20, (29) 

for a squeezed state with a squeezing angle of zero degree. Here, we define the 
squeezing level as the base 10 logarithm of squeezing (i.e. the noise variance in dB). 

The variance of a pure squeezed state, Vq, and the impurity, ri, are experimentally 
determined by the measured squeezing, Vsqz, anti-squeezing, Va-sqz, and the 
corresponding homodyne detector efficiency, t/hd [13], 

Vo = 7^^% (30) 

_ r/Hj,(2 - Vsgz - Va-sqz) - (1 - Kg^)(l - Vg^sqz) .^iN 

""'" i'2-Vsqz-Va-sqz)VHD ' ^ ' 

A typical squeezed vacuum state after the 'magic' refiector in a Schrodinger kitten 
state generation experiment as shown in figure 1 can be obtained with r2 = 0.08, 
squeezing of Vsqz = 0.661 (—1.8 dB), anti-squeezing of Va-sqz = 1-995 (-1-3 dB), and 
a homodyne detection efficiency of tjhd = 68%. According to (30) and (31), the 
corresponding pure squeezing and impurity are Vq = 0.341 (—4.67 dB) and rtotai = 
0.2438 {rtotai are the total impurity caused by ri and r2), respectively. By taking r2 into 
account, we actually have ri = 0.1771. 
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4-2. Quantum non- Gaussian character witness for an impure squeezed state and 
Schrodinger kitten state 

Based on the model developed in Section 2, a one-photon-subtracted impure squeezed 
vacuum state (Vq = -4.67 dB and impurity of ri = 0.1771 ) prepared with r2 = 0.08, 
a non-photon-number-resolving detector with Pdc = 1 x 10~^, and r]APD = 5% was 
constructed. 
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Figure 2. (a) Photon number distribution and (b) Wigner function of an impure 
squeezed vacuum state with Vb = —4.67 dB and ri = 0.1771 
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Figure 3. (a) Photon number distribution and (b) Wigner function of a projected 
state prepared with a non-photon-number-resolving detector and Vq ~ —4.67 dB, ri 
= 0.1771, ra = 0.08, Pdc^lx 10"'', tjapd = 5%, mo = 100%, and mode purity = 1. 



The photon number distribution and Wigner function of the input impure squeezed 
vacuum state and the projected state are shown in figures 2 and 3. We can see that 
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both states possess positive Wigner functions at the origin. However, the projected 
state is clearly a non-Gaussian state. Figure 4 shows the impure squeezed vacuum state 
with anti-squeezing operation. We observed that the impure squeezed vacuum state is 
a non-classical but Gaussian state, as expected. 
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Figure 4. (po, Pi) trajectory of an impure squeezed vacuum state with anti-squeezing 
operation when anti-squeezing parameter s varies from 0.18 to 0.57 (black dot-dash hue 
with circle). The dot-dash green and solid pink lines represent the classical boundary 
and the quantum Gaussian boundary, respectively. The dash blue line overlaps with 
the solid pink line and corresponds to the physical limit, po + pi = 1. 
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Figure 5. Characteristic identification of the projected state, (a) (po, Pi) trajectory 
of the projected state with anti-squeezing operation when anti-squeezing parameter 
s varies from to 1 (dot-dash black line). The dot-dash green and solid pink 
lines represent the classical and quantum Gaussian boundaries. The dash blue 
line corresponds to the physical limit, pq + pi = 1. (b) The optimal witness 
W{aopt,s) - Wciaopt) 



As indicated in figure 5(a), the one-photon projected Schrodinger kitten state 
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(corresponding to s = 0) is located beyond the classical state boundary but is within 
the quantum Gaussian state boundary. Apparently, this indicates that the state is 
non-classical and Gaussian. However, after the anti-squeezing operation, the kitten 
state crosses the quantum Gaussian state boundary, which implies the quantum non- 
Gaussian characteristic of the kitten state (i.e. it cannot be prepared by merely mixing 
Gaussian states). Moreover, the quantum non-Gaussian character witness implies that 
quantum non-Gaussian states possess strong non-classicality. Therefore, we confirm that 
the quantum non-Gaussian character witness proposed in [18] demonstrates a powerful 
ability to identify non-classical and non-Gaussian quantum states. 

4-3. Physical mechanism underpinning each experimental imperfection 

The physical mechanism underpinning experimental imperfections are summarised in 
table 1, from which we can categorise the underlying impacts as: 1) stochastic mixture 
of an m-photon-subtracted state with an M-photon-subtracted squeezed vacuum state, 
where (M > m), 2) photon number redistribution of an ideal Schrodinger kitten state, 
and 3) statistical mixture of the projected state with unprojected states. 

Table 1. Physical mechanism underpinning experimental imperfections 

Impacts Experimental imperfections Physical mechanism 

1) APD inefficiency, ?7apd Statistically mixed m-photon-subtracted 
Non-photon-number-resolving ability squeezed vacuum state with M-photon- 
Squeezing level, Vq (M > m) subtracted state (i.e. m-click 
Reflectivity, r2, events actually result from m+1, m+2 

... -photon subtraction) 

2) Input state impurity, ri N (N > 0) photon in an ideal kitten 
Homodyne detection inefficiency, r}^^, state is redistributed to N-1, N-2, 

N-3,...0 

3) APD dark count, Pdc Statistically mixed projected state with 
Mode impurity, s' unprojected state (i.e. m click events 

may actual originate from m-1, m-2 
.... -photon subtraction) 



To illustrate the above statements, figure 6 shows the photon number distributions 
for states prepared with (a) an IMPNRD with r]APD = 50%, and Pdc = 0, and (b) a 
perfect NPNRD with tjapd = 100%, and Pdc = 0, where all other related parameters 
are set to be the same. The photon number distributions are quite similar, which 
implies that the advantage of photon number resolution dramatically reduces when the 
detection efficiency is low. 

As another example, figure 7 gives the photon number distributions for Schrodinger 
kitten states prepared with (a) an input squeezed vacuum state with an impurity 
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Figure 6. Photon number distribution of Schrodingcr kitten states prepared with (a) 
an inefficient APD and (b) a non-photon-number-resolving APD . In both cases, other 
parameters are: Vq ~ —A.&I dB, ri = 0, r2 = 0.08, Pdc = 0, and mode purity — 1. 



of 0.1771 and 100% homodyne detection efficiency, and (b) an pure input squeezed 
vacuum state but 80% homodyne detection efficiency, where other related parameters 
are the same. The similarity in photon number distributions indicates the equivalent 
quantitative effect of the input state impurity and homodyne detection inefficiency on 
the projected states. Furthermore, figure 8(a) and (b) verifies that a high dark count 
probability of an APD and mode impurity demonstrate equivalent impacts on kitten 
states. 
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Figure 7. Photon number distributions of Schrodinger kitten states with (a) an impure 
input state and (b) inefficient homodyne detection. In both cases, other parameters 
are: Vq = —4.67 dB, r2 = 0.08, Pdc = 0, 'qAPD = 1, and mode purity = 1. 
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Figure 8. Plioton number distributions of Schrodinger kitten states prepared with 
(a) an APD with dark counts and (b) mode impurity. In both cases, other parameters 
are: Vo = -4.67 dB, n = 0.1771, ra = 0.08, f]APD = 1, and rjHD = 1- 



4-4- Dependency of the quantum non-Gaussian character witness and Wigner function 
at origin, W(0,0), on experimental imperfections 

Although superconducting transition edge sensors (TESs) with photon-number- 
resolving abihty are available [13, 14], commercially available APDs are still widely 
used as photon-number detectors in Schrodinger kitten state generation experiments 
since cryogenic environments are required for TESs. Typical dark count probabilities 
and detection efficiencies of commercially available Si-APDs (for 860 nm from Perkin 
Elmer Ltd.) and InGaAs-APDs (for telecommunication wavelength from ID Quantique 
Ltd.) are listed in table 2. Si-APDs perform better than InGaAs-APDS due to their 
lower dark count probabilities and higher detection efficiencies. A group of typical values 
for other related parameters are shown in table 3. We will focus on discussing the 
impacts of experimental imperfections on the quantum non-Gaussian character witness 
value and Wigner function of a one-photon-subtracted vacuum state prepared with a 
Si-APD(SPCM-AQR-12) and an InGaAs-APD (id200). 



4.4- 1- Effects of squeezing level, r2, and tjapd Figure 9 shows the variation of the 
quantum non-Gaussian character witness value and W(0,0) for a Schrodinger kitten 
state with the squeezed vacuum state level, Vq in dB, prepared with 1) perfect photon- 
number-resolving detector (PNRD), 2) a perfect non-photon-number-resolving detector 
(NPNRD), 3) an imperfect photon- number- resolving detector (IMPNRD), and 4) an 
imperfect non-photon-number- resolving detector (IMNPNRD). The other parameters 
used in these simulations to generate figure 9 are listed in table 3. Figure 9(a) shows 
that a negative Wigner function can not be observed for a Schrodinger kitten state 
at telecommunication wavelengths prepared with an imperfect photon-number detector 
and the given experimental parameters. However, it is easy to obtain a quantum non- 
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Table 2. Comparison of Si-APDs & InGaAs-APDs performance 
Detector type Dark count probability(P(ic) Quantum efficiency (r/^p^i) 



Si-APD(SPCM-AQR-12) 5x10"^ 45% 

Si-APD(SPCM-AqR-13) 2.5x10"^ 45% 

Si-APD(SPCM-AQR-14) 1 x 10"^ 45% 

Si-APD(SPCM-AQR-15) 5x10-"^ 45% 

Si-APD(SPCM-AqR-16) 2.5x10"'^ 45% 



InGaAs-APD(id200) 1 x 10"'' 10% 

InGaAs-APD(id220) 1 x 10"^ 10% 

-5 



(under different settings) 2.5 x 10"^ 15% 



5 X 10-5 20% 



Gaussian state once the squeezing level (i.e. the minimum noise variance) exceeds —0.8 
dB. On the contrary, when the squeezing level is as small as —0.4 dB, it is possible to 
obtain a Schrodinger kitten state with a negative Wigner function when prepared with 
a Si-APD , as shown in figure 9(b). 

In addition, it can be seen that there is an optimal value for the squeezing level 
of the input state to obtain maximal character witness value and minimal W(0,0). 
This is because when the squeezing level is lower than this optimal value, as shown 
in figure 9(a), not enough photons are subtracted. Consequently, the dark counts 
will dominate the 'real' click events. For large levels of squeezing (e.g. —6 dB), the 
probability to subtract more than one photon is dramatically increased. As a result, the 
one-photon- subtracted squeezed vacuum state is contaminated by two or three-photon- 
subtracted states. Therefore, the optimization of the input squeezing level is critical in 
the experiment design. 

Furthermore, a perfect PNRD demonstrates a significant advantage over a perfect 
NPNRD for both the Si-APD and InGaAs-APD cases. However, such an advantage of 
the PNRD disappears in the case of the InGaAs-APD when imperfections, such as dark 
count and detection inefficiency, are taken into account. 

As both the dark count probability and detection efficiency of the InGaAs-APD 
are inferior to those of the Si-APD, we investigated the impacts of squeezing level 



Table 3. Typical experiment parameters used in the simulation 

Parameter Typical value 

Squeezing level, Vq —4.67 dB 

Input state impurity, ri 0.1771 

Reflectivity, r2 0.08 

Mode purity, s' 0.8 

Homodyne detection efficiency, tjhd 85% 
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Figure 9. Quantum non-Gaussian character witness & W(0,0) vs input state 
squeezing level, Vb, for (a) an InGaAs-APD and (b) a Si-APD. Dash lines and solid 
lines represent 1/F(a,s)— Waia) on left vertical axis and W(0,0) on right vertical 
axis, respectively. Red: PNRD, Green: NPNRD, Pink: IMPNRD, Blue: IMNPNRD. 
Pink and blue lines overlapped in (a). PNRD is a perfect photon- number- resolving 
detector and NPNRD is a perfect non-photon-number-resolving detector. IMPNRD is 
an imperfect photon-number-resolving detector and IMNPNRD is an imperfect non- 
photon- number-resolving detector. 
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Figure 10. Squeezing level, Vq, impacts under (a) different detection efficiencies (10%, 
45% and 65% ) with P^c == 1 x 10"'' and (b) dark count probabilities (1 x 10"'', 5 x 10-^ 
and 5 x 10^^) with i]apd = 10%. Dash and solid lines represent W{a, s) — WG{a) on 
left vertical axis and W(0,0) on right vertical axis, respectively. Light blue: IMPNRD 
with (a) r^APD = 65% and (b) P^^ = 5 x 10"^, Black: IMNPNRD with {a)-qAPD = 
65% and (b) P^c = 5 x 10-^ Green: IMPNRD with (a) t^apd = 45% and (b) P^c 
= 5 x 10-^ Red: IMNPNRD with (a) t^apd = 45% and (b) P^c = 5 x 10-^ Pink: 
IMPNRD with (a) t^apd = 10% and (b) P^^ = 1 x 10"'', Blue: IMNPNRD with (a) 
r]APD = 10% and (b) P^^ = 1 x 10^''. 
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for InGaAs-APDs with different detection efficiencies and dark count probabilities 
(see figure 10). Figure 10(a) implies that the advantage of an IMPNRD becomes 
distinguishable at higher squeezing levels when the APD detection efficiency is enhanced 
to 45%. While very little difference is observed between an IMPNRD and an IMNPNRD 
in figure 10(b), despite the dark count probability being reduced to the same level as the 
Si-APD {Pdc = 5xl0~^). This reveals that the low detection efficiency of InGaAs-APDs 
substantially diminishes the advantage of an IMPNRD over an IMNPNRD. 

More importantly, reducing the dark count probability from 1 x 10~^ to 5 x 10~^ 
results in a significant decrease of W(0,0), which indicates that a lower dark count 
probability is critical to obtain negative Wigner function in Schrodinger kitten state 
generation. Comparing figure 10(a) and (b), we can see that W(0,0) obtained from a 
TES with Pdc = 1 X 10~^ and rjTES = 65% is similar to that from an InGaAs-APD 
with Pdc = 5 X 10"^ and tjapd = 10 %. This confirms the experimental result in 
reference [14], and implies that a low dark count probability is more infiuential than a 
higher detection efficiency or photon-number-resolving ability in one-photon-subtracted 
squeezed state generation experiments at telecommunication wavelengths. Therefore, 
for an InGaAs-APD with adjustable detection efficiency and dark count probability , 
such as the detector id220, the setting with lowest dark count probability is preferable 
despite the smaller detection efficiency. 
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Figure 11. Quantum non-Gaussian character witness & W(0,0) vs r2 for (a) an 
InGaAs-APD and (b) a Si-APD. Dash hues and sohd hues represent W{a, s) — Wa{a) 
on left vertical axis and W(0,0) on right vertical axis, respectively. Red: PNRD, Green: 
NPNRD, Pink: IMPNRD, Blue: IMNPNRD. 



As discussed in table 1, figure 11 shows that the refiectivity, rg, of the 'magic' 
refiector has a similar impact to the input state squeezing level, Vq, due to the same 
physical mechanism. The refiectivity, r2, must be optimized to obtain the maximum 
quantum non-Gaussian character witness value and minimum W(0,0). The optimal 
value of r2 for a Schrodinger kitten state prepared with a Si-APD (r2opt=0.01) is notably 
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Figure 12. Quantum non-Gaussian character witness & W(0,0) vs APD efficiency for 
(a) an InGaAs-APD and (b) a Si-APD. To obtain negative-valued Wigner function, 
teh requirements for the InGaAs-APD and Si-APD required to ijapd > 28%, and 
VAPD > 2.5%, respectively. Dash lines and solid lines represent W{a,s) — Wc{a) on 
left vertical axis and W(0,0) on right vertical axis, respectively. Red: PNRD, Green: 
NPNRD, Pink: IMPNRD, Blue: IMNPNRD. 



smaller than that prepared with an InGaAs-APD (r2opt=0.09). Under the circumstance 
of a Si-APD, the smaller r2 results in a larger character witness value and a deeper 
W(0,0). However, if r2 is too small, then it is easy to induce false clicks since the 
number of real APD counts, which are proportional to r2, will be lower than the amount 
of dark counts. Therefore, it is necessary to compromise a small r2 that is still higher 
enough to ensure the count rate is larger than the dark count rate. This has been 
validated by the results reported in most kitten state generation experiments using 
Si-APDs [10, 11, 12, 13]. 

The effects of the APDs quantum efficiencies are shown in figure 12(a) and (b). As 
expected, the quantum efficiency has a similar impact on the quantum non-Gaussian 
character witness value and W(0,0) to squeezing level and r2- It is noted that an 
IMPNRD does not demonstrate superiority to IMNPNRD until the detection efficiency 
increases to a specific value (for instance for the InGaAs-APD and Si-APD, r]APD = 
7.5% and 4%, respectively), which strengthens the argument that APD inefficiency 
and non-photon-number resolving ability give equivalent effects. Furthermore, when 
the detection efficiency is too low, the performance of an NPNRD is superior to an 
IMPNRD. This trend is due to the IMPNRD suffering from both inefficient detection 
and dark counts. However, when the detection efficiency is increased (for the InGaAs- 
APD and Si-APD, tjapd = 70% and 50%, respectively), the impact from dark counts 
dominates in the IMPNRD. As a result, the performance of IMPNRD is gradually 
superior to that of the NPNRD, and approaches the performance of a perfect PNRD. 
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4.4-2. Input state impurity, ri, and homodyne detection efficiency, tjhd As shown in 
figures 13 and 14, the same physical mechanism underlying the input state impurity 
and homodyne detection inefficiency in Schrodinger kitten state generation results 
in similar quantitative impacts on the quantum non-Gaussian character witness and 
W(0,0). The superior performance of the Si-APD to the InGaAs-APD culminates in 
lower requirements on the purity of the input state and homodyne detection efficiency 
at ~ 860 nm kitten generation experiments compared to such states generation at ~ 
1550 nm. 
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Figure 13. Quantum non-Gaussian character witness & W(0,0) vs input state 
impurity ri for (a) an InGaAs-APD and (b) a Si-APD. To obtain negative- valued 
Wigner functions for the InGaAs-APD and the Si-APD, ri are required to be less 
than 0.15 and 0.20, respectively. The corresponding requirements to obtain quantum 
non-Gaussian states are ri < 0.38 and ri < 0.47 for the InGaAs-APD and the Si-APD, 
respectively. Dash lines and solid lines represent W{a, s) — Wcia) on left vertical axis 
and W(0,0) on right vertical axis, respectively. Red: PNRD, Green: NPNRD, Pink: 
IMPNRD, Blue: IMNPNRD. Pinks and blue lines overlapped in (a). 



4.4-3. Dark count probability of a photon-number detector, Pdc, and mode impurity 
Both dark count probability and mode impurity cause 'false' clicks (i.e. click event is 
recorded even if no photon is actually subtracted). As shown in figure 15, to obtain 
W(0,0) < 0, the dark count probability of an InGaAs-APD is required to be less than 
2 X 10~^, which is one order of magnitude lower than that of the Si-APD {Pdc < 2 x 10"^) 
due to the lower detection efficiency of the InGaAs-APD. However, the dark count 
probabilities of most commercially available photon-number detectors for 1550 nm are 
far larger than that for 860 nm, as shown in table 3. This gives a sound reason as to why 
it is difficult to obtain negativity in the Wigner function for Schrodinger kitten states 
at telecommunication wavelengths. 

The change seen in the character witness and W(0,0) for varying mode purity is 
shown in figure 16. The PNRD is superior to the NPNRD when the mode purity is 
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Figure 14. Quantum non-Gaussian character witness & W(0,0) vs homodyne 
detection efficiency for (a) an InGaAs-APD and (b) a Si-APD. To obtain negative- 
valued Wigner functions for the InGaAs-APD and the Si-APD, tjhd are required to 
be higher than 0.88 and 0.83, respectively. The corresponding requirements to obtain 
quantum non-Gaussian states are tjhd > 0.61 and 7]hd > 0.53 for the InGaAs-APD 
and the Si-APD, respectively. Dash lines and solid lines represent W{a, s) — Wc{a) on 
left vertical axis and W(0,0) on right vertical axis, respectively. Red: PNRD, Green: 
NPNRD, Pink: IMPNRD, Blue: IMNPNRD. Pink and blue lines overlapped in (a). 




(a) InGaAs-APD 



Figure 15. Quantum non-Gaussian character witness & W(0,0) vs dark count for (a) 
an InGaAs-APD and (b) a Si-APD. To obtain quantum non-Gaussian states, Pdc are 
required to be less than 6 x 10""^ and 3 x 10"^, for the InGaAs-APD and the Si-APD, 
respectively. Dash lines and solid lines represent W{a,s) — WG{a) on left vertical axis 
and W(0,0) on right vertical axis, respectively. Red: PNRD, Green: NPNRD, Pink: 
IMPNRD, Blue: IMNPNRD. Pink and blue lines overlapped in (a). 
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high, but the advantage of PNRD gradually declines when the mode purity is too low 
to successfully projected the Schrodinger kitten state. 
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Figure 16. Quantum non-Gaussian character witness & W(0,0) vs modal purity for 
(a) an InGaAs-APD and (b) a Si-APD. To obtain negative-valued Wigner functions 
for the InGaAs-APD and the Si-APD, mode purity are required to be higher than 
0.83 and 0.78, respectively. The corresponding requirements for mode purity to obtain 
quantum non-Gaussian states are higher than 0.65 and 0.61, for the InGaAs-APD and 
the Si-APD, respectively. Dash lines and solid lines represent W{a,s) — Woia) on 
left vertical axis and W(0,0) on right vertical axis, respectively. Red: PNRD, Green: 
NPNRD, Pink: IMPNRD, Blue: IMNPNRD. Pink and blue lines overlapped in (a). 
Pink, blue and green lines are quite close in (b). 



5. Conclusions 



We quantitatively analysed the impacts of a full set of experimental imperfections 
on Schrodinger kitten state generation in terms of the quantum non-Gaussian 
character witness and Wigner function. According to the comparison between 
Schrodinger kitten states prepared with an InGaAs-APD and a Si-APD, the 
inferiority of telecommunication-wavelength photon-number detectors justifies the 
higher requirements on the optical experimental parameters to obtain negativity in 
the Wigner function. Furthermore, the lower detection efficiency of commercially 
available photon-number detectors dramatically degrades the superiority of the photon- 
number-resolving detector for one-photon projected Schrodinger kitten state generation 
at telecommunication wavelengths. The dark count probability of InGaAs-APDs is 
required to be on the order of 10~^ to obtain negative values at W(0,0). This discussion 
on the effects of various experimental parameters guides the analysis of kitten state 
generation experiments for particular wavelengths. It is clear that Schrodinger kitten 
state generation at telecommunication wavelengths presents numerous challenges but we 
can overcome these obstacles with thoughtful planning and careful experimental design. 
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